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1. Preface 

The purpose of these notes is to gather together several facts concerning the 
combinatorial aspects of diagram notation in field theory, as well as to say some- 
thing of the use of combinatorics in probability theory and stochastic processes. 
For simplicity we have restricted to Boson systems only. The motivation comes 
from the fact that, for the full potential of the diagram notation to be realized, it 
must be supported by appropriate mathematical techniques. In the general physics 
community, there are those who never use diagrams and those who never use any- 
thing else. In favour of diagrammatic notation it should be emphasized that it is a 
powerful and appealing way of picturing fundamental processes that would other- 
wise require unwieldy and uninspired mathematical notation. It is also the setting 
in which a great many problems are generated and settled. On the other, it should 
be noted the diagrams are artifacts of perturbation theory and typically describe 
unphysical situations when treated in isolation. The main problem is, however, one 
of mathematical intuition. The author, for one, has always been troubled by the 
idea that, in a series expansion of some Green function or other, some diagrams 
my have been omitted, or some have been wrongly included, or some are equivalent 
but not recognized as such, or the wrong combinatorial weights are attached, etc. 
Ones suspects that a great many physicists habour similar concerns and would like 
some independent mechanism to check up on this. 

Now, it is a fair comment to say that few people get into theoretical physics 
so that they can do combinatorics! Nevertheless, it turns out that any attempt to 
study statistical fields through perturbation theory becomes an exercise in manipu- 
lating expressions (or diagrams representing these expressions) and sooner or later 
we find ourselves asking combinatorial questions. The intention here is to gather 
together several combinatorial features arising from statistical field theory, quan- 
tum field theory and quantum probability: all are related in one way or another 
to the description of either random variables or observables through their moments 
in a given state. The connotation of the word combinatorics for most physicists 
is likely to be the tedious duty of counting out all diagrams to a particular order 
in a given perturbation problem and attaching the equally tedious combinatorial 
weight, if any. On the contrary, there is more combinatorics than mindless enumer- 
ation and the subject develops into a sophisticated and elegant process which has 
strong crossover with field theoretic techniques. It is hoped that the notes delve 
just enough into combinatorics, however, so as to cover the essentials we need in 
field theory, yet give some appreciation for its power, relevance and subtly. 

John Gough 

Nottingham 
2005 



1. PREFACE 

Frequently used symbols 

B n Bell numbers (number of partitions of a set of n items) . 

e (/) exponential vector for test function /. 

T + (f)) the Bose Fock space over a Hilbert space f). 

h n number of hierarchies on a set of n items. 

ft (X) set of all hierarchies of a set X. 

the Mobius function 
^3 (X) set of all partitions of a set X. 
^3„ set of all partitions of {1, 2, • • • n}. 
tyn,m set of all partitions of ty n consisting of m parts 
tyf (X) set of all partitions of a set X finer than {X}. 
*P C (X) set of all partitions of a set X coarser than {{x} : x e X}. 
s (n, m) Stirling numbers of the first kind. 
S (n, m) Stirling numbers of the second kind. 
& n set of permutations on {1, • • • , n}. 

&n,m set of permutations on {1, • • • , n} having exactly m cycles. 
W t the Wiener process. 

x^ falling factorial power. 

x T rising factorial power. 

© direct sum. 

£g> tensor product. 

Cg) symmetrized tensor product. 
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Combinatorics 



One of the basic problems in combinatorics related to over-counting and under- 
counting which occurs when we try to count unlabeled objects as if they where 
labeled and labeled objects as if they were unlabeled, respectively. It is curious 
that the solution to an early questions in probability theory - namely, if three 
dice were rolled, would a 11 be more likely than 12? - was answered incorrectly 
by under-counting. The wrong answer was to say that both were equally likely 
as the are six ways to get 11 (6+4+1=6+3+2=5+5+1=5+4+2=5+3+3=4+4+3) 
and six ways to get 12 (6+5+1=6+4+2=6+3+3=5+5+2=5+4+3=4+4+4). The 
correct solution was provided by Pascal and takes into account that the dice are 
distinguishable and so, for instance, 6+5+1 can occur 3!=6 ways, 5+5+1 can occur 
3 ways while 4+4+4 can occur only one way: so there are 6+6+3+6+3+3=27 
ways to get 11 and only 6+6+3+3+6+1=25 ways to get 12. The wrong solution 
for dice, however, turns out to be the correct solution for indistinguishable sub- 
atomic particles obeying Bose statistics. The same over-counting of microscopic 
configurations is at the heart of Gibbs' entropy of mixing paradox in statistical 
mechanics. 

In this section, we look at some basic enumerations and show how they are of 
use in understanding moments in probability theory. 

1. Combinatorics and Moments 

1.1. Ordinary Moments. Let X a random variable and let F (x) = Pr {X < x}. 

The expected value of some function / of the variable is then given by the Sticltjes 
integral E [/ (X)] = / (x) dW (x) . We assume that the variable has moments 
to all orders: that is E [X n ] exists for each integer n > 0. The moment generating 
function is then defined to be 

(1.1) M (t) = — E PH t n = E [e tx ] 

and so M (t) = e tx d¥ (x) is Laplace transform of the probability distribution. 
The moments arc recovered from M (t) by differentiating to the appropriate order 
and evaluating at t = 0, vis. 

(1.2) E[X n ]=^M(t) 

ai t=o 

1.2. Factorial Moments. Suppose we now work with the parameter z = e* — 
1, then M (t) = E (1 + z) x and we consider Taylor expansion about z = 0. Now 

A ,.n „u„„„ (x\ _x{x- 1) (x - 2)---(x-n+l) 

n! 



note that (1 + zf = X)^Lo i^j^ wnere {^j 



l 
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It is convenient to introduce the notion of falling factorial power 

(1.3) x ln :=x(x-l)(x-2)---(x-ra+l) 
as well as a rising factorial power 

(1.4) x T " := x (x + 1) [x + 2) . . . (x + n - 1) . 

(We also set x L ° = x T ° = 1.) 
It then follows that 

oo 

(1.5) E [z x ] = ^2 — E [X in ] z n 

n=0 

and we refer to E [Xl™] as the n-th falling factorial moment: clearly E [z x ] acts 
as their generating function. 

1.3. Stirling's Numbers. We now derive the connection between ordinary 
and falling factorial moments. We begin by remarking that the right-hand sides 
of (|1.3|) and (|1.4|) may be expanded to obtain a polynomial in x of degree n with 
integer coefficients. The Stirling numbers of the first and second kind are defined 
as the coefficients appearing in the relations 

(1.6) x Tn = ^s(n,m)x m ; x n = ^ 5 (n, m) x lm 

m m 

Evidently, the Stirling numbers of the first kind are integers satisfying s (n, m) > 
0. It turns out that this is also true of the second kind numbers. We also 
have s(n,m) = = S (n,m) for m ^ l,...,fi. It is easy to see that x^ n = 
J2 m (~ l)" +m s (n, m) x m and x n = J2 m ( — l) n+m S (n, m) x^ m and from this we 
see that the Stirling numbers are dual in the sense that 

(1.7) ]T (-If +fe a (n, k) S (k, m) = S nm . 

k 

The Stirling's numbers satisfy the recurrence relations (Stirling's Identities) 

(1.8) s(n + l,m) = s (n, m — 1) + ns (n, m) ; 

(1.9) S(n + l,m) = S (n,m — 1) + mS (n,m) 

with s(l,l) = 1 = 5(1,1). [Evidently s(n,n) = S(n,n) = S (n, 1) = 1. The 
relations x^ n+1 ^ = x^ n x(x + n) and ixi im = (x — m + m) x 1 -" 1 = i^ m+1 '+mi im 
lead to the identities when we equate coefficients.] 

This means that the Stirling numbers may then be generated recursively using 
a construction similar to Pascal's triangle, vis. 



s (n, m) S (n, m) 



n\m 
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n\m 


1 


2 


3 


4 


5 


6 
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120 


274 


225 
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90 
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From relation we see that the ordinary and falling factorial moments are 
related by 

E [X n ] = £ (n, m) E [X lrn ] , 

m 

(1.10) E[X ln ] = ^(-l)" +m s(n,m)E[X m ]. 

m 

1.4. Cumulant Moments. Cumulant moments K n are defined through the 
relation YlnLi — K nt n = InM (t) or 



^-E[X«]r = exp 

n=0 ' Ln=l ' J 



and one sees from expanding and comparing coefficients of t that 

E [X] = K%, 

E[X 2 ] = K 2 + l4, 

E [X 3 ] = n 3 + 3k 2 «i + «i, 
E [X 4 ] = k 4 + 4k 3 ki + 3^2 + 6k 2 Ki + k 4 , 
etc., 

or inversely 

ki = E [X] , 

k 2 = E [X 2 ] - E [X] 2 

n 3 = E [X 3 ] - 3E [X 2 ] E [X] + 2E [X] 3 

Hi = E [X 4 ] - 4E [X 3 ] E [X] - 3E [X 2 ] 2 + 12E [X 2 ] E [X] 2 — 6E [X] 4 , 
etc. 

1.5. Examples. 

Example 1 (Standard Gaussian). We take F (x) = (27r) 1 ^ 2 e~ y2 / 2 dy lead- 
ing to the moment generating function 

M{t) = e t2 ' 2 . 

We see that all cumulant moments vanish except k% = 1. Expanding the moment 
generating function yields 

( (2k)\ 

{ 0, n = 2k + l. 



Example 2 (Poisson). We take X to be discrete with Pr{X — n} — ^jA"e~ A 
for n — 0, 1,2, The parameter A must be positive. The moment generating 

function is readily computed and we obtain 

M{t) = exp{A(e*-l)}. 

Taking the logarithm shows that all cumulant moments of the Poisson distribution 
are equal to A: 

K n = A, for all n = 0, 1, 2, • • • . 
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Likewise, M (t) = exp {Xz} where z = e* — 1 and so we see that the falling factorial 
moments are 

E [X ln ] = A". 

The ordinary moments are more involved. From <|i._?fl|) we see that they are poly- 
nomials of degree n in A : 

(1.11) E[X n ]=J2S(n,m)X m 

m 

and in particular A is the mean. We may expand the moment generating function 
through the following series of steps 

oo 1 

M (t) = £-V( e '-i) m 



m! 

m— 



^m~k ^kt 



m=0 fc=0 v 7 

oo ^ m / \ oo ^ 



m! 1 — ' \ k / ' — ' n 

m=0 k=0 x 7 ra=0 



and comparison with ll j.J j|) yields the Stirling identity 

*<».»> -si; 



mi 

k=0 



~ k k n . 



Example 3 (Gamma). Let X &e a positive continuous variable with ¥ (x) = 
r (A) Jq y x ~ 1 e~ y dy , (the Gamma function is T (A) = J °° y x ~ 1 e~ y dy), the its 
moment generating function is 

M(t) = {i-ty x . 

Now (1 — £) A = X)^Lo f ^ ( — = S^Lo TIT^™^™ so it s momen t s are 
(1.12) E[r] = A l "=^s(n,m)r. 

2. Fundamental Enumerations 

We now want to draw attention to the fact that the various families of numbers, 
appearing in the formulas for the moments of our main examples above, have an 
importance in combinatorics |15| . 

(2k)\ 

For instance, the factor —. — - occurring in the moments of the Gaussian are 
well-known as the count of the number of ways to partition a set of 2k items into 
pairs. The Stirling numbers also arise as counts of classes of permutations and 
partitions, as we shall see below. 

2.1. Pair Partitions. A pair partition of the set {1,2,..., 2k} consists of k 
pairs taken from the set so that every element of the set is paired with another. 
We shall denote the collection of all such pair partitions by V^k- Evidently, \V2k \ = 
(2k)\ 

2ki\ ' we nave 2fc x (2k — 1) choices for the first pair, then (2k — 2) x (2k — 3) 
for the second, etc. This gives a total of {2k)\ however we have over-counted by a 
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factor of k\, as we do not label the pairs, and by 2 k , as we do not label the elements 
within each pair either. It is convenient to set \V2k+i \ = since clearly we cannot 
partition up an odd number of elements into pairs. 

2.2. Permutations. The set of permutations, 6„, over {1, ...,n} forms a 
non-abelian group under composition. We shall use the notation er° = id, a 1 — 
a, a 2 = a o a, etc. 

Given a permutation a £ 6 n , the orbit of a number i £ {1, . . . , n} under a is the 
sequence \i, a (i) , a 2 (i) , . . . } . As the orbit must lie within {1, . . . , n} it is clear that 
a k (i) = i for some k > 0: the smallest such value is called the period of the orbit 
and clearly orbit repeats itself beyond this point (a n+k (i) — <j n (i)) . The ordered 
collection \i, a (i) ; a 2 (i) ; . . . ; a fc_1 (i)] is referred to as a cycle or more explicitly a 
k- cycle. Cycles will be considered to be equivalent under cyclic permutation in the 
sense that [x±; X2', ■ ■ ■ ; Xk] is not distinguished from [x2j X3; . . . ; Xk',xi], etc. Thus 
each fc-cycle is equivalent to k sequences depending on where on the orbit we choose 
to start. Clearly orbits arising from the same permutation a either coincide or are 
completely disjoint; this simple observation leads to the cyclic factorization theorem 
for permutations: each permutation a can be uniquely written as a collection of 
disjoint cycles. 

Lemma 1. Let & n ,m be the set of permutations in & n having exactly m cycles. 
Then the number of permutations in & n ,m is given by the Stirling numbers of the 
first kind 

(2.1) \6 n , m \ = s(n,m). 

Proof. This is proved by showing that |© n , m | satisfies the same recurrence re- 
lation (|1.8I) as the Stirling numbers of the first kind, that is |© n +i,ro| = |©n,m-i| + 
n \& n ,m\ -Now |© n +i, m | counts the number of permutations of {1, 2, . . . , n + 1} hav- 
ing m cycles. Of these, some will have n + 1 as a fixed point which here means 
that [n + 1] is unit-cycle: as we have m — 1 cycles left to be made up from the 
remaining labels {1, 2, . . . , n}, we see that there are |S„ im _i| such permutations. 
Otherwise, the label n + 1 lies within a cycle of period two or more: now if we 
take any permutation in & n .m then we could insert the label n + 1 before any one 
of the labels i G {1, . . . , n} in the cyclic decomposition - there are n |6 njm | such 
possibilities and the second situation. Clearly |©i,i| = 1 = |6 n ,n| while |6„, m | = 
if m > n. Therefore |6„. m | = s (n, m). | 

2.3. Partitions. Let X be a set. We denote by *}3 (X) the collection of all 
partitions of X, that is, A = {Ai, . . . , A m } £ *|3 (X) if the Aj are mutually-disjoint 
non-empty subsets of X having X as their union. The subsets Aj making up a 
partition are called parts. The set X is trivially a partition - the partition of X 
consisting of just one part, namely X itself. All other partitions are called proper 
partitions. 

If A" = {1, . . . ,n} then the collection of partitions of X will be denoted as 
while the collection of partitions of X having exactly m parts will be denoted 

as yftn.m, • 
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Lemma 2. The number of partitions of {1, . . . 
Stirling number of the second kind S (n, m): 

(2.2) |«p n , m | =S(n,m) 



n} having exactly m parts is the 



Proof. To prove this, we first of all show that we have the formula |^P n +i,m| = 
|^Pn,m-i| + m I^Pn,m|- This is relatively straightforward. We see that |?JJ n +i, m | 
counts the number of partitions of a set X = {1, . . . , n, n + 1} having m parts. 
Some of these will have the singleton {n + 1} as a part: there will be |^P n , m _i| of 
these as we have to partition the remaining elements {1, . . . ,n} into m — 1 parts. 
The others will have n + 1 appearing with at least some other elements in a part: 
we have |^p n ,m| partitions of {1, . . . , n} into m parts and we then may place n + 1 
into any one of these m parts yielding m|^3„ !m | possibilities. Clearly = 1 

and |$JJn,n| — 1 while |*P„. m | = if m > n. The numbers |^ n ,m| therefore satisfy 
the same generating relations (|1.9() as the S (n, m) and so are one and the same. | 

As a corollary, we get the following result. 

The total number of partitions that can be made from n symbols, termed the 
n-th Bell number and denoted by B n , is given by 



(2.3) B n = |?JJ n | 

The first few Bell numbers are 



m— 1 



S (n, m) 



n 


1 


2 


3 


4 


5 


6 


7 


8 


B„ 


1 


2 


5 


15 


52 


203 


877 


4140 



For instance, the set {1,2,3,4} can be partitioned into 2 parts in 5(4,2) = 7 
ways, vis. 



{{1,2}, {3, 4}}, 
{{1,3}, {2, 4}}, 
{{1,4}, {2, 3}}, 

and into 3 parts in S (4, 3) = 6, vis. 

{{1},{2},{3,4}} : 
{{1},{4},{2,3}}. 
{{2}, {4}, {2, 4}}. 



{{1},{2,3,4}}, 
{{2}, {1,3, 4}}, 
{{3}, {1,2, 4}}, 
{{4}, {1,2, 3}}, 

{{1},{3},{2,4}}, 
{{2}, {3}, {1,4}}, 
{{3}, {4}, {1,2}}. 



2.3.1. Occupation Numbers for Partitions. Given A € ?fi(X), we let rij (A) 
denote the number of parts in A having size j. We shall refer to the rij as occupation 
numbers and we introduce the functions 



N(A) = J2 n {A) ,E{A) = Y, jnj (A) . 



If A" = {1, . . . ,n} then the collection of partitions of X will be denoted as 
^P„. while the collection of partitions having exactly m parts will be denoted as 

Note that ifA£ <P„, m then N (A) = m and E (A) = n. 

It is sometimes convenient to replace sums over partitions with sums over 
occupation numbers. Recall that a partition A will have occupation numbers 
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n = (m, ri2, n.3, . . . ) and we have N (A) = m + ni + n.3 + • • • and i5 (-4) = 
rti + 2ri2 + 3n 3 + • • • . 

We will need to count the number of partitions with E (A) = n leading to the 
same set of occupation numbers n; this is given by 

1 77' 

(2.4) p(n 



! . . . 



(2!)^ (3!)" 3 • m!n 2 !n 3 ! 

The argument is that there are n\ ways to distribute the n objects however we 
do not distinguish the nj parts of size j nor their contents. We remark that the 
multinomials defined by 



{Z1,Z 2 , ■■■)= P( U ) Z " l2: 2 S 



are known as the Bell polynomials. 

2.3.2. Coarse Graining and Mobius Inversion. A partial ordering of *p (X) is 
given by saying that A =4 B if every part of A is a union of one or more parts of B. 
In such situations we say that A is coarser than B, or equivalently that B is finer 
than A. The partition consisting of just the set X itself is coarser then every other 
partition of X . Likewise, the partition consisting of only singletons is the finest 
possible. 

Whenever A =4 B we denote by rij (A\B) the count of the number of parts of A 
that occur as the union of exactly j parts of B. We also introduce the factor 

M (^i/3)=-n{-^- i ) ! }" 3( " 4l6) - 

3>l 



Theorem 1. Let : ^P(<-f) 1— > C be given and let a function be defined by 
$ (B) — X^^t^B ^ (-4)- ^ e relation may be inverted to give 

Proof. Essentially we must show that 

£ v(A\B) = S A , c . 

A4B^C 

Suppose that Aj £ A, then we may write Aj as the union of kj, say, parts of B 
and rj, say, parts of C. Evidently we will have 1 < rj < kj. By considering all the 
possible partitions of these kj parts of C (for each j = 1, I) we end up with all 
the partitions B finer than C but coarser than A. The sum above then becomes 

N(A) ( k 3 

II E^ 1 ^ to -1)15 

however, observing that (r — 1)! = s (r, 1) and using the duality of the first and 
second kind Stirling numbers l|1.7|l . we see that this is proportional to J^J - =1 5i,kj- 
This gives the result. | 
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Note that if ip is a function of the subsets of X and if 4> (A) = J^bca V* (-^0 
then we have the relation ip (A) — ^2 B(ZA (— l)'" 4 '"' 5 ' 0(B) which is the so-called 
Mdbius inversion formula |13| . The above result is therefore the corresponding 
Mdbius inversion formula for functions of partitions. 

2.4. Hierarchies. Let X be a set and let 

tyf (X) : = "set of all finer partitions of X" , 
ty c (X) : = "set of all coarser partitions of X" 

where we mean (X) excluding the coarsest, finest partition respectively. 

A hierarchy on Xis a directed tree having subsets of X as nodes, where A is 
further along a branch from B if and only if A c B, and where X is the root of the 
tree and {x} ,i6l, are the leaves (terminal nodes). 

For instance, let X = {1, • • ■ ,6} then a hierarchy is given by taking three nodes 
A = {2, 3}, B = {1, 2, 3} and C = {4, 5, 6}. 

X 




{1} {2} {3} {4} {5} {6} 

There are two equivalent ways to describes hierarchies, both of which are useful. 
2.4.1. Bottom-up description. Let us consider the following sequence of parti- 
tions: 

AM = (AfU'U'^fW: 4 1J ={1},4 1) = {2,3}, 4 1} = {4,5,6}; 

-4 (2) - {4 2) ,4 2) } g r : 4 2) = {^U^} , 4 2) = {4 13 } ; 
= {4 3) }er(-4 (2) ) : 4 3) = {4 2) ,4 2) }. 

This equivalently describes or example above. 

In general, every hierarchy is equivalent to such a sequence 

A {1) e r (*) , A (2) e r , A (3) e r („4 (2) ) , • • • 

and as the partitions are increasing in coarseness (by combining parts of predeces- 
sors) the sequence must terminate. That is, there will be an m, which we refer to 
as the order of the hierarchy, such that A^ = j^™ -1 )}. Evidently, m measures 
the number of edges along the longest branch of the tree. 
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In our example, we can represent the hierarchy as 

H = 4 3) 

= {4 2 \4 2) } 

= {{4 1) ,4 1) },{4 1) }} 

= {{{1}, {2, 3}}, {{4, 5, 6}}}. 

This is an order three partition - each of the original elements of X sits inside three 
braces. 

2.4.2. Top-down description. Alternatively, we obtain the same hierarchy by 
first partitioning {l,--- ,6} as «8W = jsf^B^j where b{ 1} = {1,2,3} and 

B^ = {4,5,6}, then partitioning B^ as {1} and {2,3}, and finally partitioning 
all the parts at this stage into singletons. 

In general, every hierarchy can be viewed as a progression: 

Q3 (1) G % f (X) , G W f (B) : B G <8 (1) } , {«8jf G V f (-B) : B G Q5 (2) } 

Eventually, this progression must bottom out as we can only subdivide X into finer 
partitions so many times. Again the maximal number of subdivisions is again given 
by the order of the hierarchy. 

We shall now introduce some notation. Let Sj (X) denote the collection of all 
hierarchies on a set A. We would like to know the values of h n , the number of 
hierarchies on a set of n elements. We may work out the lowest enumerations: 

When n = 2 we have the one tree /\and so h 2 — 1. When n = 3 we have the 
topologically distinct trees 




and, when we count the number of ways to attach the leaves, we have h 3 = 1 + 3 = 4 
possibilities. 

When n = 4 when have the topologically distinct trees 




which implies that h A = 1 + 4 + 6 + 12 + 3 = 26. We find that 



n 


1 


2 


3 


4 


5 


6 


7 


8 


K 


1 


1 


4 


26 


236 


2752 


39208 


660032 



This sequence is well known in combinatorics and appears as sequence A000311 on 
the ATT classification. It is known that the exponential generating series h (x) = 
J2 n ^h n x n converges and satisfies exph (x) = 2h{x) — x + 1. 

3. Prime Decompositions 

There results of this section are really necessary for our discussion on field 
theory, but have been included for completeness. As is well known, the primes 
are the indivisible natural numbers and every natural number can be decomposed 
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uniquely into a product of primes. Each natural number determines an unique 
sequence of occupation numbers consisting of the number of time as particular 
prime divides into that number. In our discussion below, we meet a calculation 
of functions that can be defined by their action on primes. Here we encounter 
an argument for replacing a sum (over occupation sequences) of products with a 
product of sums (over individual occupation numbers). This argument will recur 
elsewhere. 

It turns out that the cumulant moments play a role similar to primes insofar 
as they are indivisible, in a sense to be made explicit, and every ordinary moment 
can be uniquely decomposed into cumulant moments. 

3.1. The Prime Numbers. The natural numbers are positive integers N = 
{1, 2, 3, • • • }. If a natural number m goes into another natural number n with no 
remainder, then we say that m divides n and write this as m\n. Given a natural 
number n, we shall denote the number of divisors of n by d(n), that is 

d (n) := # {m : m\n} . 

Likewise, the sum of the divisors of n is denoted as 

k\n 

k 

A natural number, p, is prime if is has no divisors other than itself and 1. That 
is, d{p) = 2. The collection of primes will be denoted as V and we list them, in 
increasing order, as 

Pi = 2, p 2 = 3, p 3 = 5, pi = 7, p 5 = 11, etc. 



Theorem 2 (Prime Decomposition Theorem). Any natural number m can be 
uniquely decomposed as a product of primes: 

oo 

(3-1) m = H( Pj r. 

The numbers nj = nj (m) give the number of times the j th prime, pj, divides 
an integer m. In this way, we see that there is a one-to-one correspondence between 
the natural numbers and the collection of "occupation numbers" n — {nj)°° =l where 

have < Y^jLi n j < °°- 

Theorem 3 (Euclid). There are infinitely many primes. 

PROOF. Suppose that we new the first N primes where N < oo, we then 
construct the number q = YijLi Pj + 1- If we t f y t° divide q by any of the known 
primes pi, ■ ■ -Pn, we get a remainder of one each time. Since any potential divisor 
of q must be factorizable as a product of the known primes, we conclude that q has 
no divisors other that itself and one and is therefore prime itself. Therefore the list 
of prime numbers is endless. | 
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If to and n have no common factors then we say that they are relatively prime 
and write this as to _L n. The Euler phi function, (p (n) , counts the number of 
natural numbers less than and relatively prime to n: 

if (n) := # {to : to < n, to _L n} . 



3.2. Dirichlet Generating Functions. Let a — (a n )?L , be a sequence of 
real numbers. Its Dirichlet generating function is defined to be 



oo 



n— 1 

Let a and 6 be sequences, then their Dirichlet convolution is the sequence c = a*dC 
defined by 

m|n 

^ ^ Q>mbn/ va- 
ra 

Lemma 3. Let a and b be sequences, then Da^b = ®c where the sequence c is 
the Dirichlet convolution of a and b. 

Proof. 

oo oo , oo oo , oo 

a b ~2^2^ n s- m s~2^2^ ( nrn ) s ~ 2.^ k s ' 

n—1 ra—1 n—1 m—1 k— 1 



Let / be a function defined on the natural numbers. The function is said to be 
multiplicative if 

/ (nm) = / (n) f (to) , whenever mln. 

If we furthermore have / (nm) — f(n)f(m), for every pair of natural numbers, 
then the function is said to be strongly multiplicative. For instance, / (n) = n s is 
strongly multiplicative. However, the Euler phi function is multiplicative, but not 
strongly so. 

Lemma 4. Let f be a multiplicative function of the natural numbers. Then 

pEV ra>Q F 

Proof. This is a consequence of the unique prime decomposition l|3.1|) for 
any natural number. If / is multiplicative and to = rXj=i (Pj)™ 3 ' then / (to) = 
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Il^li / (Pj J )- The Dirichlet generating function for the sequence / (to) is then 
^ / \ f (to) 



TO" 

m>l 



n j=l 

oo oo 



j=l n=0 

- nz- 

pEV n>0 y 



The replacement ^2 n YiJLi 9 ( n j) ^ T\°r^i S^Lo 9 ( n ) use d above is an elemen- 
tary trick which we shall call the J2 Y\ <-> Y\ 2 trick. (It's the one that is used to 
compute the grand canonical partition function for the free Bose gas!) 

Lemma 5. Let f be a strongly multiplicative function of the natural numbers. 
Then 



»/<•>- n(i-^)' 



Proof. If / is strongly multiplicative, then / (p n ) = f (p) n . We then encounter 
the geometric sequence 

v ./•(?") _ v /(p)" _ ( 1 f(p)Y 1 

' i n ns / j n ns \ 



n>0 P ^0 P V P 



3.3. The Riemann Zeta Function. The Riemann zeta function, ((s), is 
the Dirichlet generating function for the constant sequence 1 = (l)^L r That is, 



oo 

(3-2) C(s):=£-. 



n=l 



We note that ( (0) and C (1) are clearly divergent. However, £ (2) = Y^=i ^ 
— . It is clear that ( (s) is in fact an analytic function of s for Re (s) > 1. 

Lemma 6 (Eulcr). For Rc (s) > 1, 

(3.3) c( S )=nf 1 -^ 



pev 



i 

p" 



Proof. This follows from the observation that ( (s) = T>i (s) and that the 
constant sequence 1 corresponds to the trivially strongly multiplicative function 
f(n) = l. I 
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An immediate corollary to the previous lemma is that £ 2 (s) — X>d (s) = 
d{n) 

2-n=l n s ■ 

3.4. The Mobius Function. The Mobius function is the multiplicative func- 
tion n determined by 

/*(p n ) ■- 
for each prime p. 

Lemma 7. £> M = ^. 

proof. = e„>i^? - n P£P L>o^ = n P ev(i-p- s ) = 



1, 


n = 0; 


-1, 


n = 1; 


o, 


n > 2. 



Lemma 8. Lef (6n)^Li fee a gwen sequence and set a n = Y^k ™ ^fc> 



fe 

PROOF. We evidently have that £> a = £> h C and so D fc = C _1 £>a = SV^V | 

As an application, we take b n — n. Then a„ = X)fe ™ ^ which is just the sum, 
s (n), of the divisors of n. We deduce that n = Ylk A 4 (f ) s 



CHAPTER 2 



Boson Fock Space 



Why do the Gaussian and Poissonian distributions emerge from basic enumer- 
ations? We shall try and answer this by looking at Bosonic field and their quantum 
expectations. 

1. Boson Statistics 

Identical sub-atomic particles arc indistinguishable to the extent that there 
states must be invariant under arbitrary exchange (or, more, generally permutation) 
of their labels we might attach to them. Suppose wc have a system of n identical 
particles, then it happens in Nature that one of two possibilities can arise: either 
the total wave function ip (xi,x 2 , • • • ,x n ) is completely symmetric, or completely 
anti-symmetric under interchange of particle labels. The former type are called 
bosons and we have 

V> (x<y(l),Xa(2), ' ' ' ,X a (n)) = V> ( x l , x 2 , ■ ■ ■ , x n ) , for all (7 G 6„, 

while the latter species are called fermions and we have 

tp {xa(i),x a{ 2), ■ ■ ■ ,x a{n) ) = (-i f ip(x 1 ,x 2 , ■■■ ,x n ), for all a e S„, 

where (— l) a denotes the sign of the permutation. (Here, the coordinates Xj give 
all necessary labels we might attach to a particle: position, spin, color, etc.) 

1.1. Fock Space. Let () be a fixed Hilbert space and denote by f)®" the closed 
linear span of symmetrized vectors of the type 

(1.1) Vl® • • • ®<Pn ■= ^ Yl ^(1) ® • • • ^(n)- 

' itE6„ 

where ipj G t). To understand better the form of f)® n , let {ej : j = 1, • • • } be an 
orthonormal basis of f). Suppose that we have n particles with n\ particles are in 
state ei, n 2 in state e 2 , etc., then there is accordingly only one total state describing 
this: setting n = {rij : j = 1, 2, • • • }, the appropriately normalized state is described 
by the vector -known as a number vector - given by 




n\ factors ri2 factors 

The span of such states, with the obvious restriction that . rij = n, yields f)®". 

The Fock space over () is then the space spanned by all vectors |n) when only 
the restriction J2j n j < 00 applies. We also include a no-particle state called the 
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Fock vacuum and which we denote by fi. The Fock space is then the direct sum 

(i.3) r + (ij) = ©r. 

n=0 

The vacuum space I)® is taken to be spanned by the single vector fi. 

Often, the number state vectors are not the most convenient to use and we now 
give an alternative class. The exponential vector map e : f) i— > T + (f)) is defined by 

(1-4) s(f) = ®^ (4=^ n 



with /®" the n-fold tensor product of / with itself. The Fock vacuum is, in partic- 
ular, given by fi = e (0). We note that (e (/) \e (g)) = J2n>a h. = CX P (/l^), 
whence the name exponential vectors. (These vectors are called Bargmann vector 
states in the physics literature, while their normalized versions are known as coher- 
ent state vectors.) The set e (f)), consisting of all exponential vectors generated by 
the test functions in (), is linearly independent in T + (h). Moreover, they have the 
property that e (S) will be dense in the Boson Fock space whenever S is dense in 
f). 

1.2. Creation, Annihilation and Conservation. Boson creation and an- 
nihilation fields are then defined on r + (fj) by the actions: 

B + {4>) fi®---®f n : = Vn+1 (f>®fi®- ■ ■ <8>/ n ; 

B~{4>) fi® ■••<£>/„ : = -)= ^ (<l>\fj) fi® • • • ®fj® • • • ® f n - 

They satisfy the canonical commutation relations 

[s-($,fl + (V0] = (<^>, 

(1.5) [B-(0),B-(^] = 0= [B+(<f>),B+M]. 

Likewise, let M be an operator acting on (). Its differential second quantization 
is the operator dT + (M) acting on Y + (f)) and defined by 

dT+ (M) (fi® ■ ■ ■ ®f n ) := (M/i) ® ■ ■ ■ ®/„ + ■ ■ ■ + /i® ■ ■ ■ ® (M/„) . 

2. States 

2.1. Gaussian Fields. Consider (fi| 5 £ « (/„) • • • B £ ^ (f 2 ) B e « (/i) fi), where 
the e's stand for + or — , this is a vacuum expectation and we may compute by 
the following scheme: every time we encounter an expression B~ (fi) B + (fj) we 
replace it with B + (fj)B~ (fi) + (fi\fj)- The term (fi\fj) is scalar and can be 
brought outside the expectation leaving a product of two less fields to average. Ul- 
timately we must pair up every creator with an annihilator otherwise we get a zero. 
Therefore only the even moments are non-zero and we obtain 
(2.1) 

n 

E {n\B^(hn)---B^(f 2 )B^(fi)Sl)= ]T II 
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Here fe, <7j)j=i ^ T~ > 2n is a pair partition: the pj correspond to annihilators and 
the qj to creators so we must have pj > qj for each j; the ordering of the pairs is 
unimportant so for definiteness we take q n > ■ ■ ■ > > q\. We may picture this 
as follows: for each i E {1, 2, . . . , 2n} we have a vertex; with B + (fi) we associate a 
creator vertex with weight fi and with B~ (fi) we associate an annihilator vertex 
with weight fi. A matched creation-annihilation pair (jpj, qj) is called a contraction 
over creator vertex qj and annihilator vertex pj and corresponds to a multiplicative 
factor (f Pj \f qj ) and is shown pictorially as a single line, 



annihilator vertex at p 



creator vertex at q 



contraction between vertices p and q 



9 



We then consider a sum over all possible diagrams describing. 
Setting all the fj equal to a fixed test function / and let Q (/) = B + (f)+B~ (/) 
then we obtain 

n| q (f) k n) = (n\ [b+ (/) + b~ (/)] k n) = \\f\\ k \v k \ . 



The observable Q (/) therefore has a mean-zero Gaussian distribution in the Fock 
vacuum state: 



For instance, we have \V±\ — = 3 and the three pair partitions {(4,3) (2, 1)}, 
{(4, 2) (3, 1)} and {(4, 1) (3, 2)} are pictured below. 

r\r\ QQ rTPi 

432 1 432 1 432 1 



We remark that (|2.1|l is the basic Wick's theorem in quantum field theory |5], 
and can be realized in terms of Hafnians j2j. However, the result as it applies to 
multinomial moments of Gaussian variables goes back to Isserlis jllj in 1918. 



2.2. Poissonian Fields. More generally, we shall consider fields B ± (.) on 
some Fock space on r + (h). Consider the vacuum average of an expression of the 
type 

(fi| B+ (f n ) a(n) B- (g n f in) ...B+ (/i) Q(1) B~ { 9l ) m Vi) 

where the a/s and (3j's are powers taking the values or 1. This time, in the 
diagrammatic description, we have n vertices with each vertex being one of four 
possible types: 
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scattcrer emitter absorber constant 




The typical situation is depicted below: 



1 . ^7i?K7f^vh 



i(5) i(4) i(3) z(2) i(l) 

Evidently we must again join up all creation and annihilation operators into 
pairs; we however get creation, multiple scattering and annihilation as the rule; 
otherwise we have a stand-alone constant term at a vertex. In the figure, we can 
think of a particle being created at vertex i (1) then scattered at i (2) , i (3) , i (4) 
successively before being annihilated at i (5). (This component has been highlighted 
using thick lines.) Now the argument: each such component corresponds to a unique 
part, here {i (5) , % (4) , i (3) , i (2) , i (1)}, having two or more elements; singletons 
may also occur and these are just the constant term vertices. Therefore every such 
diagram corresponds uniquely to a partition of {1, . . . , n}. Once this link is made, 
it is easy to see that 

(n| B + (f n ) ain) B~ (g n f {n) ..-B+ (A) Q(1) B~ ( 9l ) m 

q,/3g{o,i}" 

(2.2) = ^ II (fli(fc)l/i(fc-l)>---(ffi(3)|/i(2)>(fli(2)|/i(l)>- 

AeVr, {i(k)>->i(2)>i(i)}eA 
If we now take all the fj and gj equal to a fixed / then we arrive at 

<0| (/) + !) (S- (/) + !)]" O) = E E H/ll""™ 

m=ore<p„ >m 



m) ||/|| 2( "- ro) 



m=0 



Note that a part of size fc contributes ||/|| 2 ^ fc ^ so a partition in ty n ,m with parts 
of size k u ...k m contributes ||/|| 2^+-+ k ™- m ) = ||/|| 2( "" m) . 
It therefore follows that the observable 

(2-3) N{f):=(B + (J-]+\\f\\\(B-' 1 



has a Poisson distribution of intensity ||/|| 2 in the Fock vacuum state. 
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2.3. Exponentially Distributed Fields. Is there a similar interpretation 
for Stirling numbers of the first kind as well? Here we should be dealing with cycles 
within permutations rather than parts in a partition. Consider the representation 
of a cycle (i (1) , i (2) , . . . , i (6)) below: 




To make the sense of the cycle clear, we are forced to use arrows and therefore we 
have two types of lines. In any such diagrammatic representation of a permutation 
we will encounter five types of vertex: 




♦ 



An uncontracted (constant) vertex indicates a fixed point for the permutation. 

Let us consider the one-dimensional case first. The above suggests that we 
should use two independent (that is, commuting) Bose variables, say bf,bf. Let 
us set 

&±/:= 6*012 + 1!®^ 
then we see that b + will commute with b~ . We note that 

b + b~ = b^b^ (8)12 + 12® &2"&2 + &r ® &2 + b i ® bt + 1 
and here we see the five vertex terms we need. 

Let ft i and ft 2 be the vacuum state for and respectively, then let ri = £7i(g) 
f?2 be the joint vacuum state. We wish to show that N has a Gamma distribution 
of unit power (an Exponential distribution!) in this state. 

First of all, let Q = b + + b~ so that Q = Qi <g> 1 2 + li Qi where Q 3 = b+ +bj. 
Now each Qj has a standard Gaussian distribution in the corresponding vacuum 

state e tQ i Clj) = e* 2 / 2 )and so 

(Q\ e tQ ft) = <Qi| e tQl fti) (n 2 | e tQ2 Q 2 ) = e { ' \ 

Therefore <fi| Q 2 " ft) = 2" |P 2n | = However Q 2n = ^ m O (b + ) m (b-) 2n ~ m 

(remember that b+ and b~ commute!) and so (ft| Q 2n ft) = ( 2 ™) (ft| (b+) n (b~) n ft) 
which should equal ( 2 ™) (ft| (6+6")" ft). Therefore 

(ft| (&+&-)" ft) = n! 
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and so, for t < 1, 

(n| exp{tb+b-} n) = 
Therefore N = b + b~ has an exponential distribution in the joint vacuum state. 

The generalization of this result to Bosonic fields over a Hilbert space t) is 
straightforward enough. First we need the notion of a conjugation map on Hilbert 
spaces: this is a conjugate-linear map j : f) h- > h. For instance, let {e„} n be a fixed 
complete basis for f) then an example of a conjugation is given by j (^2 n c n e n ) = 

c n e n- 



Theorem 4. Let I) be a separable Hilbert space and j an conjugation on [). Let 
B (.) be operator fields define on T + (h) ® T + (f)) by 

(2.4) B± (/) = S± (/) ® 1 2 + li ® Bf (j7) 

where the B i (.) are usual Bosonic fields on the factors T + (h), i (= 1,2). T/ien 
[B~ (/) , £+ (#)] = for all f,g€h and if N (/, g) := B+ (/) B~ (g) i/ien we have 
the following expectations in the joint Fock vacuum state Q = Oi £g> f^2 : 

(/«,$») •••JV"C/i,si) n) 

( 2 -5) = X] II (ffi(fc)l/i(fe-l)) • • * (fl'*(2)l/t(l)> X <ffi(l)|/i(fc)> 

<re6„ (i(l),i(2),...,i(fe))e<7 

The proof should be obvious from our discussions above so we omit it. The 
sum is over all permutations er G & n and each permutation is decomposed into 
its cycles; the product is then over all cycles (i (1) , i (2) , . . . ,i (fc)) making up a 
particular permutation. We note that the representation corresponds to a type of 
infinite dimensional limit to the double- Fock representation for thermal states . 

2.4. Thermal Fields. There is a well known trick for representing thermal 
states of Bose systems. To begin with, let f) = L? (R 3 ) be the space of momentum- 
representation wave functions, then the thermal state ( • )p is characterized as the 
Gaussian (sometimes referred to as quasi- free), mean zero state with 



(B+(f)B-(f))^ = J \f(k)\ 2 g^(k)d 3 k 

where 

(We have the physical interpretation of (3 as inverse temperature, fj, as chemi- 
cal potential, and E (k) the energy spectrum function.) Let us denote by Qp 
the operation of pointwise multiplication by the function on h, that is, 

( W) = Bf># (*) f We check that ( t B+ (/) + B ~ (/)] 2 )„ = /) 

Pi ft 

with Cfj tll = 2g fj41 + 1 = coth - ^ ~ ^ . 

Theorem 5. Let I) be a separable Hilbert space and j an conjugation on fj. Let 
B ± (.) be operator fields define on T + (fj) <X) T + (f)) by 



(2.6) B {f) = Bf y/g+lf ) ® l a + li ® S 2 T 0' V?/) 
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where the B { (.) are usual Bosonic fields on the factors T+ (()), i (= 1, 2), and g is a 
positive operator on fj . Then the fields satisfy the canonical commutation relations 
and their moments in the joint Fock vacuum state CI = Ox ® ^2 <we precisely 
the same as for the thermal state when we take g — Qp 

To see this, note that 

(n| e 4^ + (/)+^-(/)] n) = (o x | e <K(/)+Brw] 0l ) (n 3 | e *K(/)+^-(/)] n 2 ) 
= exp j-i (/| ( e + 1) /)| cxp |-i (/| ( e + 1) /)| 

= exp|-i(/|C/)| 

where C = 2g + 1 is the covariance of the state. The first factor describes the 
so-called spontaneous emissions and absorptions while the second factor describes 
the so-called stimulated emissions and absorptions. In the zero temperature limit, 
P — > oo, we find that g — > and so the second factor fields (stimulated emissions 
and absorptions) become negligible. 
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Field Theory 

1. Introduction 

Let us suppose that we have a fixed space $ of functions over a set A. A function 
tp = {tp x : x e A} in $ will be called a field realization. By a field <f) we mean an 
observable taking values in $ (that is to say, the field realization are "eigen-values" 
for the field) and we suppose that we are supplied with state which we shall denote 
as ( • ). Formally, we may think of the field as a family of operators 4> — {<fi x : x e A} 
and we suppose that \_<j> x , 4> y ] = for all x,y <G A. Since we are assuming that the 
field is commutative, we may in fact revert to classical probability and think of the 
field as being a random variable taking values in The expectation of a (typically 
nonlinear) functional F = F [</>] of the field is to be understood as an integral over 

(1.1) TO]>= / F [tp] dv (tp) , 

where v is the probability measure over $ corresponding to the state. We shall 
frequently use the notation ( • ) = when we want to emphasize the correspon- 
dence. 

The label x € A is assumed to give all relevant information such as position, 
spin, etc. When A is a finite set, then the mathematical treatment is straightfor- 
ward. Otherwise, we find ourselves having to resort to infinite dimensional analysis. 
We also introduce a dual space 3 of fields J = {J x : x E A} which we call the source 
fields. Our convention will be that the source fields carry a 'contravariant' index 
while the field, and its realizations, carry a 'covariant' index. The duality between 
fields and sources will be written as 

In the case where A is a finite set, say of cardinality N, a realization tp <E $ can 
be viewed as a set of N numbers and so we can identify <i>, and likewise 3, as M N . 
In this case tp x ,J x means X^xeA fx^ x an< ^ so the notation just implies an Einstein 
summation convention over A. If dtp x denotes standard Lebesgue measure, and let 
S [•] be some functional, called the action, such that S = J $ exp {S [tp]} dtp < oo, 
then, we have the finite-dimensional probability measure, v, on $ determined by 

dv(<p) = ±e s M ]Jd Vx . 
" xeA 

The general situation which we are really interested is where A is continuous. If 
we want A = M. d , then we should take Z to be the space of Schwartz functions on R d 
and $ to be the tempered distributions. (The field realizations being more singular 
than the sources!) Here, the duality is denoted is (tp, J) = J Rd tp x J x dx and we again 
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shorten to just <p x J x . We shall therefore adopt an Einstein summation/integration 
convention over all indices from now on. The appropriate way to consider random- 
izing the field in the infinite-dimensional case, will then be to consider probability 
measures on the Borel sets of Z, that is, on the c-algebra generated by the weak 
topology of the Schwartz functions 1 . 

2. Field Calculus 

In this section, I want to present a way of handling functional derivatives ex- 
ploiting the commutativity of multiple point derivatives. The notation that results 
has similarities to one introduced by Guichardet [8] for investigating Bosonic Fock 

spaces. We shall denote by the functional derivative wrt. J x . For X a finite 

oJ x 

subset of A we shall adopt the notations 

II m = n^- 

x£X xGX 

(In both cases, we need only the set X - the ordering of the elements is irrelevant!) 
We of course have -j-j- (J v ) — S x . Similarly, we shall write ip x for Yl xe \<p x , an d 

sm n s 



6<p x Lm Sif x 

2.1. Analytic Functionals. A multi-linear map T : x™3 i— > C is called a 
tensor of covariant rank n and it will be determined by the components T Xl ... Xn such 
that T (J(i), ■ • • , J( n )) = T Xl ... Xn J x ^ • ■ • Likewise, we refer to a multilinear 

map from x™$ to the complex numbers as a tensor of contravariant rank n. 

A functional F = F [ J] is said to be analytic in J if it admits a series expansion 
of the form F [J] = J2n>o hfxi—x n J Xl ■ • • J x " where f Xl —x„ are the components 
of a completely symmetric covariant tensor and as usual the repeated dummy in- 
dices are summed/integrated over. A more compact notation is to write the series 
expansion as 

x 

$ $ 

It is easy to see that j-^ ■ ■ ■ j-p^F [J] = ~£ n >o 7J.fyi-y m x 1 -x n J Xl • • ■ J Xn and 
this now reads as 




2.2. The Leibniz rule. Functional derivatives obey all the standard rules of 
calculus, including the Leibniz rule and we have the natural extension to multiple 
derivatives: 

S \X\ _ S \Y\p S \X/Y\q 

(2 ' 1} jjx( FG )-2-, x ^jv- ^jxpr- 



At this stage we shall stop with the functional analysis details and proceed as if every thing 
is well-defined. The mathematically inclined reader can fill in the details, while everyone else 
could well live without them. 
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The generalization of this for several factors is 

X\X\ > S \Yi\ p S \Y m \p 

SjX {*l---*rn)-2^ Yi ,.. Ym&x SjYl SJ Y m 

where the sum is over all partitions of X into m labelled parts (as opposed to the 
unlabeled ones we have considered up to now). 

2.3. Differential Formulas. We now derive a useful result. 
Lemma 9. Let W = W [J] then 

Proof. This is easily seen by induction on n. As — — e w = e w ^^-, the 
identity is true for n = 1. Now assume that it is true for n, then 

v n slAl]V w — v J n slAlw 

SJ X 8J x6 e SJ X ^ M <W A ^ SJ X M ^ 

.4e<p(x)Ae_4 _4e<p(Jf) Ae.4 

however, the first term on the right hand side is a sum over all parts of X U {a;} 
having x occurring as a singleton, while the second term, when differentiated wrt. 
J x , will be a sum over all parts of X U {x} having x in some part containing at 
least one element of X. Thus we may write the above as 

r5l*l +1 . „ X\A\ W 

_£ P w — w sr ft £ !£_ 

«5jXU{x} e 2^ 11 W A • 

-4e<P(XU{a;}) Ae.4 

The identity the follows by induction. | 

2.4. Legendre Transforms. Suppose that W = W [ J] is a convex analytic 
function. That is, W is a real-valued analytic functional with the property that 

(2.3) W [tJi + (1 - 1) J 2 ] < tW [Ji] +(l-t)W [J a ] 

for all < £ < 1 and Ji,J% £ 3- The Legendre-Frenchel transform of TF is then 
defined by 

(2.4) I>]= inf {IF[J]-(<^, J)}. 

r [ip] will then be a concave (i.e. — T is convex) analytic functional in tp and we may 
invert the formula as follows: 

(2.5) W [J\ = sap {T[<p] + (<p, J)}. 

If the functional W is taken to be strictly convex, that is, if we have strict 
inequality in (|2.3() . then the infimum is attained at a unique source J — J [<p] 
for each fixed ip, and so T [ip] = W [J[ip]\ — (ip, J [ip]). Moreover, we may invert 
J : $ i— > Z to get a mapping : -3 i— > <f>, and for fixed J the supremum is given by 
<j) [J] and so IF [J] = T \6 [J]\ + \ 4> [J] , J)- The extremal conditions are then 

« = _^T_ - _ <W 
" 6<p m > f'-SJ'- 
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Let W" [J] be the symmetric tensor with entries W'^ y [J] 



S 2 W 



This will 



SJ x 5Jy ' 

be positive definite - it will be interpreted below as the covariance of the field in 
the presence of the source J. Likewise, if we let T" [<p] be the linear operator with 

entries — , then we nave 



r" xy \u>] = 



SJy' 
SJ y 



and so we conclude that W" [J] and —T" [J] = V" \4> [J]] will be inverses for each 
other. In other words, 



(2.6) 



S 2 W S 2 T 



SJ x SJy 5ip y 8<p z 



= -51 



<p=4>\A 



Lemma 10. Let F : $ R be a functional and suppose F : Z 
given by F [J] := F [cj> (J)] then 



is then 



S 2 W SF 



SF 



SJ X SJ x SJy Sip 



Proof. This is just the chain rule, as - — = 



S 2 W 



SJ X SJ x 5Jy ' 



Lemma 11. Let us introduce the tensor coefficients 

6 \x\ T 



(2.7) r [J] := 



6ip x 

then they satisfy the differential equations 



4>[J] 



AeV(Y) \AeA / 

where each za is a dummy variable associated with each component part A and 
Z_4 = {za ■ A £ ^l}. There is, as usual, an implied contraction over the repeated 
S lAl+1 W 

Za's in T XuZa and the -r— -, — r . 

5J Aij { z a} 

Proof. We shall prove this inductively. First note that for F = T x , we have 
— ^— T x = - — ^— f Xu i z } by the previous lemma. Assuming that the relation (|2.8|) 

oJy 5J y SJ z 
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is true for a set Y, let y be a free index then 



Ajy Z^ 11 



sjvsjy sjy A ^ (Y) 

Z^ | I 11 £jAu{* A } I ^ <5J{f,*} I 11 SJ^*-*} 

Ae<p(Y) I VAgA / \AeA 

Z^ I 11 §jAU{z A } I 

A6<P(YU{i/}) \AeA / 

and so the relation holds for Y U {y}. | 



Theorem 6. We have the following recurrence relation 

(29) £!5^= v ^ fMuz, f n 

1 j <5J* U M ~ Z. SJV6JP I 11 s J*u{*a} ) ' 

Aeqsf(x) \AeA / 

where, again, each za is a dummy variable associated with each component part A 
and Za = {za ■ A £ A}. 

Proof. Taking - — ^ of (|2.6|l and using the multi-derivative form of the Leibniz 
6 J 

rule, we find 

5J X \SJ*6Jy J 

Z^ ^jru{x,i/} Z^ I 11 jjau{ za } J ' 

rex ^eqj(x/Y) \aea / 

Now the Y" = X term in the last summation yields - — =r—. — T T^ y ' z ' , which will 

5J xu i x >y} 

be the highest order derivative appearing in the expression, and we take it over to 

S 2 W 

the left hand side, we then multiply both sides by — — which is the inverse of 

oJ x oJ z 

p{y,«} ) finally we get the expression 

S lXl+2 W = v v 5 2 W 6™ +2 W ( „ S^W \ f{p , q}uZA 
$jxu{x,y} Z^ Z^ SJySJP 5J Yu i x >i} \ J-L ,5JAU{^} I 

YCX A£<P(X/Y) \AeA I 

We note that the sum over all yuji} (where Y C X, but not Y — X) and partitions 
of X/Y can be reconsidered as a sum over all partitions of X U {a;}, excepting the 
coarsest one. Let us do this and set X' = X U {x} and Z' A = Za U {q}; dropping 
the primes then yields (|2.9|) . | 
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3. Green Functions 

3.1. Field Moments. The Green functions are the field moments (fi Xl • • • fi x ) , 
where X = {xi, • • • , x n } is a subset of A, and they form the components of a sym- 
metric tensor of covariant rank n = \X\ under the conventions above. We shall 
assume that field moments to all orders exist. It is convenient to introduce the 
shorthand notation 

(3.1) (<^>=(n^y 

\xeX I 

The moment generating function is given by the Laplace transform 

(3.2) Z [J] = f e^' J) dv (ip) 

J<s> 

In particular Z [0] = 1 and we have the expansion 

n>0 ' X ' '' 

The field moments are then recovered from Z [J] using functional differentiation 
S n Z 

= {4> X1 ■ ■ ■ <f> x ), or, in shorthand, this reads as 

.7=0 

. S^Z[J] 

(fix) = 



SJ x i ■■■SJ X 



6J X 



.7=0 



3.2. Cumulant Field Moments. The cumulant field moments (((fix)) (also 
known as Ursell functions in statistical mechanics) are defined through 

(3-3) (fi x ) = fl Ua)) 

Ae%5(X) A&A 

and using Mobius inversion, we get 

(3-4) «<^» = J2 »W II <<^> 

AG<$(X) a&a 

where fJ-(A) = — Ylj >1 {— (j — l)!}"^" 4 '. The first couple of cumulant moments 
are 

«*„» = <*.>, 

((<f>x<t>y)) = (fixfiy) ~ (fix) {<t>y) . 
{{<t>x<t>y<i>z)) = (fixfiyfiz) ~ (fixfiy) (fiz) ~ (<t>y<t>z) (fix) ~ (fizfix) (fiy) 
+2<0 X ><^><&>. 

Note that ((4> x )) is just the mean field while (((j> x (/) y )) is the covariance of <p x and 
4> y in the state. 

Theorem 7. The cumulant Green functions are generated through the 

(3-5) W[J] = J2-]-((<j> x )) J x =\nZ[J] 

x ' 
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Proof. For convenience, we write {CJ)j for 

= E Urn 



J*i . . . j x 3 then 



n>0 A£<Pn. AeA 

Now suppose that we have a partition consisting of rij parts of size j then the 
occupation sequence is n — (ni, 77,3, • • • ). The number of parts in the partition 
is therefore Y] j nj while the number of indices being partition is n = Y^j j n j ■ Recall 
that the number of different partitions leading to the same occupation sequence n 
is given by p (n) in equation Ij2.4|l and so 

Z W = £ ni ln 2 ln 3 \ ■ ■ • (l!)" 1 (2!)" 2 (3!)" 3 • • • {cJ) ' 2 



)7 



En n 1 

n j>\ 
IIE^T 

Il eX p{ 

7>1 ^ 



3^ 

expj] 



r 

i! 



and so Z = cxp VK. 



Note that we have used the Yl Yi Yi S trick from our section on prime 
decompositions. We could have alternatively used the formula lj2.2J) to derive the 
same result, however, the above proof suggests that we should think of cumulant 
moments as somehow being the 'primes' from which the ordinary moments are 
calculated. 

3.3. Presence of Sources. Let J be a source field. Given a probability mea- 
sure we may introduce a modified probability measure v J ' , absolutely continuous 
wrt. v 1 and having Radon-Nikodym derivative 

dv J 1 

^ M = z^j] exp{( ^ JU - 

The corresponding state is referred to as the state modified by the presence of a 
source field J. 

Evidently, we just recover the reference measure v when we put J = 0. The 
Laplace transform of the modified state will be Z v .i [K] = (eP^ ) and it is 
readily seen that this reduces to 

Z[J + K] 

Z " J [K] ~ z[j] ■ 

In particular, the cumulants are obtained through W v j [K] = W v [J + K] — W v [J] 
5 m W J [K] 

and we find {{<j> x )) v J = j^x 

K=0 
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It is, however, considerably simpler to treat J as a free parameter and just 
consider v as being the family {v J : J e j}. In these terms, we have 

(3-6) Ux)) v j = ^ [J] = E ((<W)) " JY > 

We point out that it is likewise more convenient to write 



and again we can drop the superscripts J. 

The mean field in the presence of the source, <j) [J] € $, is defined to be <f> x [J] = 
(4> x } u j and is given by the expression 

(3.7) ^ [J] = E ^ ' ' ' 0>„ JX1 ' ' ' JX " 

n>0 

and, of course, reduces to {4> x ) v when J = 0. 

3.4. States. Our basic example of a state is a Gaussian state. We also show 
how we might perturb one state to get another. 

3.4.1. Gaussian States. Let L be a linear, symmetric operator on $ with well- 
defined inverse G. We shall write g xy for the components of L and g^y for the 
components of G. That is, the equation Lip = J, or g xy if y = J x will have unique 
solution ip — GJ, or ip x — g xy J y ■ As G is positive definite, symmetric it can be 
used as a metric. It can also be used to construct a Gaussian state. 

We construct an Gaussian state explicitly in the finite dimensional case where 
|A| = N < oo by setting 

d l (</?) = , 1 = CX P \-\9 xv VxV v \ I| dip x 
J(2nf detG L J xeA 



which we may say is determined from the a quadratic action given by S 1 



\g xy <p x <p v - The moment generating function is then given by 



(3.8) Z^[J]=exp^-g xy J x J y 

In the infinite dimensional case, we may use (|3.8|l as the definition of the measure. 

The measure is completely characterized by the fact that the only non- vanishing 
cumulant is ((^i^)^ = g xy and if we now use (|3.3|l to construct the Green's func- 
tions we see that all odd moments vanish while 

( 3 - 9 ) {^xil) ■ ■ ■ ( t ) x(2k)) = ^2 9x( Pl )x( qi ) ■ ■ ■ 9x(p k )x(q k ) 

7 V 2k 

where the sum is over all pair partitions of {1, • • ■ , 2k}. The right-hand side will of 

■ , (2A:)! 
course consist of , , , terms. 

2 k kl 
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3.4.2. Perturbations of a State. Suppose we are given a probability measure /i . 



and suppose that Si [•] is some analytic functional on <I>, say Si [ip] = 
ip yi ■ ■ ■ (p Vn , or more compactly 

(3.io) s 'M = £t-U x <^. 



n>0 n\ ' 



_L,,yi---Vn 



X 



A probability measure /x, absolutely continuous wrt. /u , is then prescribed by 
taking its Radon-Nikodym to be 

^- (<p) = ^ exp {Si [ip]} 

provided, of course, that the normalization S = (exp {5/ [</?]}) < oo. 

The generating functional for n will then be 



or, more directly, 



Si[4>] e (<t>,J) 



(3.11) 



Z^, [J] = — exp <^ Si 



1 






exp 




' S ' 


{« 









S_ 

SJ 



Z, [J]- 



e (4>,J) 



We remark that, for finite dimensional fields, we may think of /j, being deter- 
mined by the action S [ip] — So [<p>] + Si [ip] where So is the action of fi . 

3.5. The Dyson-Schwinger Equation. We now derive a functional differ- 
ential equation for the generating function. 



Lemma 12. The Gaussian generating functional Z 1 satisfies the differential 
equations 

+ J X \Z 1 [J] = 



where F x [ip] = - — SL [ip]. 

°<Px 



Proof. Explicitly, wc have Z y = exp {| J x g xy J y } so that jj^Z-y — g xy J v Z 1 
which can be rearranged as 

+ J x )Z~, = 0. 



SJv 

However, S* 7 [<p] = -\g xy p x p> y and so F x [ip] = -g xy p x . | 

Lemma 13. Suppose that Z Mo satisfies a differential equation of the form 

J x \Z, o {J}=0 





' S ' 








+ J*} 


\ Mo 





and that fi, is absolutely continuous wrt. fi , is given by (<p) — A exp{S7 [</?]}. 
Then Z^ satisfies 



_5_~ 

SJ 



+ F, 



SJ 



+ J X }^[J] = 







:S2 
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where Ff 



Proof. We observe that from Q3.11|l we have 



J X Z^[J] = ~J*exp<j ,S, 
and using the commutation identity 



we find 

J x Z fi [J] 



_5_ 

SJ 



Z,„ [J] 





' 6 ' 






' S ' 




' 6 ' 




J x ,exp js 1 / 




I 


= -Ff 




exp "1 Si 




} 



— exp Sj 



_S_ 

SJ 



rz, [j\ -Ff 



' S ' 




' 8 ' 






exp "1 Si 




} z, 



-F, 



Mo 



_6_ 

SJ 



Z» [J] Ff 



S_ 
Tj 



ZAJ] 



which gives the result. | 

Putting these two lemmas together we obtain the following result 



Theorem 8 (Dyson-Schwinger) . The generating functional for a probability 
measure absolutely continuous wrt. 7 satisfies the differential equation 

' 5 ' 



(3.12) 

where F 



SS\ 



= -\g xy y y + F? 



J x \Z [J] = 



Corollary 1. Under the conditions of the above theorem, if Si [<p] = Ylx Tk\> yX Vx > 



then the ordinary moments of [i satisfy the algebraic equations 
(3-13) (^ {l)u x) = E 9xx' Ux/W}) + 9xy E W» viV}UY 



x'EX 



J XUY/ 



Proof. We now have F x = —g xy <p y +Y^x i — n v <Px- The Dyson-Schwinger 

|X|! 

equation then becomes 



+ V — v^ uY +J x Z-0 

Xjy \Y\\ V X.TY +J U 



sjy \y\\ 



6J Y 



s \x\ 



and we consider applying the further differentiation ~^jy ^° obtain 



-^(^ } ux) j+ y:^ {x}uy 

The result follows from setting J = 0. | 



s \x\ 

bxuY)j + Jjx( JXZ )= - 
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We may write (|3.13|1 in index notation as 

m 

(^x&xx ■ ■ ■ <f>x m ) = ^9xx t (4> Xl ---4>x z ---4> Xm ) 



i=l 

n>0 

where the hat indicates an omission. This hierarchy of equations for the Green's 
functions is equivalent to the Dyson-Schwinger equation. 

We remark that the first term on the right hand side of i|3.13|) contains the mo- 
ments (4>x/{x'}^j which are of order two smaller than the left hand side (^{xyux 
The second term on the right hand side of l|3.13[) contains the moments of higher 
order and so we generally cannot use this equation recursively. In the Gaussian 

case, however, we have (4>{ x }ux^j = J2 X 'ex 9xx' (4>x/{x'}^ from which we can 
deduce (|3.9[) by just knowing that (<f> x ) — and {<j> x <f> y } — g xy . 

The Dyson-Schwinger equations may alternatively stated for W, they are 

9 6Jy + ^\x\\ v 2^ H 5J A + u ' 

X 1 1 A£*p(X)A£A 

3.6. Field-Source Relations. Recall that <j> x [J] is the mean field value in 
the presence of a source J € {J and is given by 

SW 

(3.14) = ^. 

Let us assume again that the map J <—> <fi from Z to <J> is invertible and denote the 
inverse by J. That is, <f> x \J[(p]~\ = ip x . Moreover, we suppose that J admits an 
analytic expansion of the type J x [<p] = — ^2 n>0 ^T xyi '" Vn p yi ■ ■ ■ <p , or 

J x i^ = -J2wu T{x}ux ^- 

X ' 

We introduce the functional 

n>0 ' X 1 

so that 

(3.15) J* M= _il. 

°<Px 

We then recognize T [ip] as the Legendre-Frenchel transform of W [J] , that is, 
W [ J] = sup vg$ {r [ip] + {if, J)} with the supremum attained at the mean-field 

<P = 4>[J}. 

For the Gaussian state we have W 7 [J] = \g xy J x J v and we see that <j> x = g xy J v . 
The inverse map is therefore J x [ip] — g xy p y and so we obtain 

r 7 M = -\a xy v x fy 
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More generally, we have r = T 7 + Ti where 

ra>3 

where it is customary to set T xy — —g xy + n xy . Without loss of generality we shall 
assume that the lead term T x is equal to zero: this means that the state is centred 
so that 4> [0] = 0. It follows that 

(3.16) r m = (0°y - * xv ) v y - \T xvz v y Vz 

and (substituting tp — <j>) we may rearrange this to get 

(3.17) 4> x =9*y(j V +* VZ 4>z + \T yz ^ 

We note that 4> appears on the right-hand side of (|3.17fl in a generally nonlinear 
manner: let us rewrite this as <f> = GJ + f up) where / satisfies / (0) = 0. We may 
re-iterate (J3.17JI to get a so-called tree expansion 

^ = GJ + f(GJ + f(GJ + f(GJ + --- ))) 

and we know that this expansion should be re-summed to give the series expansion 
in terms of J as in l|3.7|l . 

Now T" x y [ip] = -g x y + T" xy [<p] so we conclude that 

(3.18) W" [J] ' 



L T» [J] 

where (f£ [J]) xy = Tt xy + E|x|>i jm r{x ' v}UX ^X [ J l This relation may be alter- 
natively written as the series 

(3.19) W" [J] = G 1 _ f ] J]G =G + GT" [J] G + GT" [J] Gf" [ J] G + ■ ■ ■ 

In particular, when we set J = 0, W" [0] is the covariance matrix while f " [0] = tt, 
provided we assume that the state is centred (0 = when J = 0), and we obtain 
the series expansion 

(3.20) W" [0] = G + GttG + GttGttG +■■■ 



We now wish to determine a formula relating the cumulant moments in the 
presence of the source J to the tensor coefficients T x [J]. 

Theorem 9. We have the following recurrence relation 

(3.2i) ((^uw)L = £ ((*At ftp) ^M fn((^u W } 

A£<$f(X) \AeA 

where, again, each za is a dummy variable associated with each component part A 
and Z_a — {za ■ A £ ^4}. 
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This, of course, just follows straight from ()2.9(l . The crucial thing about (|3.21|1 
is that the right hand side contains lower order cumulants and so can be used 
recursively. Let us iterate once : 



f{ P }uz A ^ f yi 



x > f^u] n w B ) 



E 



What happens is that each part A of the first partition gets properly divided up into 
sub-parts and we continue until eventually break down X into its singleton parts. 
However, this is just a top-down description of a hierarchy on X. Proceeding in 
this way we should obtain a sum over all hierarchies of X. At the root of the tree, 
we have the factor {{(fry^p)) u j an d for each node/part A, appearing anywhere in the 
sequence, labelled by dummy index za say, we will break it into a proper partition 
B e %t(A) and obtain a multiplicative factor (U^qyj) ; f^ uZe [J] where Z B 

will be the set of labels for each part B G B, 

If we set X = {xi,x%} then we have only one hierarchy to sum over and we 

find = iMX' ((KKh; ((KKK* fqziZ2 [4 

It is useful to use the bottom-up description to give the general result. First 
we introduce some new coefficients defined by 

with the exceptional case T xy := \yl>x4 l y/) j- Then we find that 
^ UX )) V J = E T {«}u2 A(m) 

>< n n ^t-i- n ^ 

(3.22) 

For instance, we have the following expansions for the lowest cumulants: 

((<f>y<i ) x 1 <i>x 2 <l>x 3 )) V J = ^yx 1 x 2 x 3 + (T„ r 'T r , ., . + •••), 
{( < Py4 > x 1 C t'x2 C f'x 3 C t'x i )) V J = r ^yx 1 x 2 x 3 x 4 + (T yXl T rX2X3X4 + •■■) + (T J/2 , l2 , 2 T,. X32 ; 4 + 

' V yxi rx2 L qx 3 Xi i J ' \ y rxix 2 ± qx 3 X4 i 

etc. 

The terms in round brackets involve permutations of the Xj indices leading to 
distinct terms. Thus, there are 1 + ^ (2) = 4 terms making up the right-hand side 
for the fourth order cumulant: The first term in round brackets corresponds to the 
hierarchy {{{^1}} , {{^2} , {^3}}} and there are 3 such second order hierarchies. 
There are 1 + (3) + 5 (1) (o) = 26 terms making up the right-hand side for the fifth 
order cumulant. 
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4. Wick Ordering 

4.1. Wick Cumulants. Up to now, we have been considering cumulant mo- 
ments of the field. Let i*i,i*2j • • • , F n be some functionals of our field, then we may 
define the object ({[Ft] ■ ■ ■ [F n ])) through 

(([Fi]---iF n ]}= e n ((n^ 

where ([F,] ■ ■ ■ [F n ]) = (F x ■ ■ ■ F n ). Thus (([F] [G])> = (FG) - (F) (G), etc. 

As the notation hopefully suggests, we treat each component in a square brack- 
ets as an indivisible whole object. For instance, we have that (([4> x ] is a 
second order cumulant (even though it has three field factors!) and it equates to 



4.2. Wick Monomials. The following definition of Wick monomial is due 
to Djah, Gottschalk and Ouerdiane j4|. It is more general than the traditional 
definition, which uses the vacuum state, and we review it in terms of our field 
calculus notation: the theorems are otherwise after |4j directly. 

Definition 1. For each finite subset X, the Wick ordered monomial : <p x : is 
the operator defined through the property that 

(4-1) (:0 X : F) = /(n [4>*\[F] 

Wxex 

for all appropriate F. 



) 



Let us first remark that the appropriate functionals F = F [(f)] are those such 
that F — F [ip] is square-integrable wrt. the measure v. Secondly, it should be em- 
phasized that the definition of : 4> x ■ depends on the choice of state or, equivalcntly, 
on v. 



Lemma 14. The ordinary field operator products 4> x can be put together from 
the Wick monomials according to the formula 



e EfnV>>) 

- z <B(x)AeA \BeA J 



vx - WX/ ' 

Ae^(x) 

Proof. First of all, we observe that we can use the short-hand form (\<Pxu{a} 

f° r ((EL ex [<Px\ [F])) by setting <j) a = F where a is an exceptional label, and we 

augment our configuration space to AU {a}. Now the ordinary moment (<Pxu{a}) 

can be expanded as S^eq3(xu{a}) IlAe^l ((^a)) an( i ^ ne partitions of X U {a} can 
be set out into two types: those that contain {a} as a singleton part, and those 
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that don't. This leads to 



■a,m/ = ( e nw)((o+ e E(n«W UQ 

\.4e$p(x) Ae^l / .4e<p(x)Ae.A \BeA / 

= <0X> + E E II ^ : *a> • 

.4e<p(x) AeA \beA I 

However, as <p a was arbitrary, we obtain the required identity. | 



COROLLARY 2. The Wick monomials satisfy the following recursion relation 
: 4>x ~ <t>x ~ (4>x) - E E f ft : ^ : • 

Theorem 10. Let X, Yi, • • • , Y m 6e disjoint finite subsets, then 
<: o V; : • • • : <^ m : <p x ) = £ [] 

AeV'(Yu--- ,Y m ;X) A £-4 

where *P (Yi, • ■ ■ , Y m ; X) is the set of all partitions o/ Yl U • • • U Y TO U X having no 
subset of any of the Yj 's as a part. 

Proof. Let n = J2jLi the result will be established by strong induction. 
The case n = 0, to begin with, is just the expansion of {4> x ) in terms of its connected 
Green's function. If the results is assumed to hold up to value n then 

(: <hi : ■ ■ ■ : </Y m+ i : <f>x) = 4>Yi : ' ' ' : 4>Y m '■ 4>Y m ux) 

- (: (f> Yl : • • • : <j) Ym : <j) x ) (^ Ym+1 ) 

(B^A 

e En 

•AeSP'OWO AeA \BeA 

and by induction up to order n we may rearrange this as 

e n«<^»- e e n«**»n«« 

Ae^' (Y lt --- ,Y m ;Y m+ iUX) AeA AeV' (Y u --- ,Y m ;X)BeVp(Y m+1 ) AeA BeB 

B=jLA 

^■ 2 ) e e e n «^» n «^c7» ■ 

Ae<pf(Y m+1 ) AeA C e<9'(Y u --- ,Y m ,A-,x) Be A cec 

The claim is that this equals X)^e<p'(Yi ••• Y m Y m i-x) TIaza ((^a))- Now the 
first term appearing in 14.211 is the sum over <P (Y%, ■ ■ ■ ,Y m ;Y m +i U X) which is 
strictly larger than *p (Yi, • • • ,Y m , Y m+1 ; X) by means of including unwanted par- 
titions of Yi U ■ ■ • U Y m+ i U X in which subsets of Y m+ i appear as parts. The second 
term of (|4.2|) subtracts all such contributions from partitions for which every part 
B that contains an element of Y m+ i only contains elements of Y m+ i. The third 
term subtracts the contribution from the remaining unwanted partitions, namely 
those where the offending parts B make up only some subset Y m+ \/A of Y m+ \. 

This establishes the result by induction. | 




Yi • • 9Y m ■■ <PA ■ 9X, 



CHAPTER 4 



Stochastic Integrals 



Stochastic processes are intermediatory between random variables and quantum 
fields. That's certainly not the view that most probabilists would have of their 
subject but for our purposes it is true enough. Indeed, there exists a natural theory 
of quantum stochastic calculus extending the classical theory of Ito to the quantum 
domain. 

We can treat stochastic processes as a special case of quantum fields over one 
time dimension. More exactly, we may view stochastic processes as regularized forms 
of white noise processes - these are stochastic processes only in some distributional 
sense, but can be realized as quantum fields theory . Instead, we shall work with the 
more traditional viewpoint used by probabilists. We first recall some basic ideas. 

We say that a random variable is second order if it has finite first and second 
moment and that a second order variable X is the mean-square limit of a sequence 
of second order variables {X n ) n if E (X — X n ) 2 — > 0. A stochastic process is a 

family of random variables {X t : t > 0} labeled by time parameter t. The process 
is second order if E [Xt] , E [A t A s ] < oo, for all times t, s > 0. 

A process {X t : t > 0} is said to have independent increments if X t — X s and 
Xt> — X s i are independent whenever the time intervals (s,t) and (s',t') do not 
overlap. The increments are stationary if the probability distribution of each X t — 
X 8 depends only on t — s, where we assume t > s. The most important examples of 
second order processes are the Wiener process {W t : t > 0} and the Poisson process 
{N t : t > 0}. The Wiener process is characterized by having increments W t — W s , 
(t > s), distributed according to a Gaussian law of mean zero and variance t—s. The 
Poisson process is characterized by having increments N t — N s , (t > s), distributed 
according to a Poisson law of intensity t — s. 

In general, given a pair of stochastic processes {X t } and {Y t }, we can try and 
give meaning to their stochastic integral X t dY t . What we do is to divide the 
interval [a, b] up as a = t < t\ < ■ ■ ■ < t n = 6 and consider the finite sum 
Y^j=o Xtj iXtj+i — Y tj ) and consider the limit in which max{tj + i — tj} — ► 0. If 
the sequence of finite sums has a mean-square limit then we call it the Ito integral 
f a XtdY t . 

There is a problem, however. The Wiener process has increments that do 
not behave in the way we expect infinitesimals to behave. Let At > and set 
AW t — Wt+At — W t , then, far from being negligible, (AW t ) 2 is a random variable 
of mean value At. The consequence is that the usual rules of calculus will not 
apply to stochastic integrals wrt. the Wiener process. Remarkably, if we work with 
suitable process there is a self-consistent theory, the Ito calculus, which extends the 
usual notions. In particular, the integration by parts formula is replaced by the Ito 
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formula, which states that 

d (X t Y t ) = (dX t ) Y t + X t (dY t ) + (dX t ) (dY t ) 

and is to be understood under the integral sign. We now present some combinatorial 
results on multiple Ito integrals: more details can be found in the paper of Rota 
and Wallstrom IK 



1. Multiple Stochastic Integrals 

Let X t be stochastic processes for j = 1, . . . , n. We use the following natural 
(notational) conventions: 

xi n > • ■ • xP = / dX^ ... I dxP = i dX^ n) ■ ■ ■ dX^ 



J[o,t] J[o,t] J[o,t] n 

In the following we denote by A™ (t) the n— simplex in [0,t} n determined by a 
permutation a £ 6„: that is, 

K (*) = {(tn, • ■ ■ , tl) € (0, t) n : t a{n) > > • • • > t a(1) } . 

We denote by A™ (t) the simplex determined by the identity permutation: that is 
t>t n > t n -i > ■ ■ ■ > ti > 0. Clearly U CTe e„A™ (t) is [0,t] n with the absence of the 
hypersurfaces (diagonals) of dimension less than n corresponding to the tj 's being 
equal. Moreover, the A™ (t) are distinct for different u. 

We also define the off-diagonal integral tl -f" to be the expression with all the 
diagonal terms subtracted out. Explicitly 

/ <">..-<<>:= £ / dX^-.-dX^. 

Take si, S2, ■ ■ ■ to be real variables and let ^4 = {^4i, . . . , A m } be a partition of 
{1, . . . , n} then, for each i 6 {1, . . . , n}, define s_4 («) to be the variable Sj where i 
lies in the part Aj . 



Lemma 15. The multiple stochastic integral can be decomposed as 

dx^ ■ ■ ■ dxP = V / dx {n \ , • • • dx (1) 



Note that there n = 1 case is immediate as the / and ^ integrals coincide. In 
the situation n = 2 we have by the Ito formula 



Jo Jo 

dXgdXg* + f dX%>dXg> + f dXWdXV 

ti>0 Jt>t 1 >t 2 >o Jo 



dX^dX^ + l dx^dxP 



i [o,t] 2 J[o,t] 
and this is the required relation. 
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The higher order terms are computed through repeated applications of the Ito 
formula. An inductive proof is arrived at along the following lines. Let x[ n+V> be 
another quantum stochastic integral, then the Ito formula is 

X (n+l)y t = I dX l n+ ^dY tn + I dX^+VdY, 
J[0,t] 2 " + J[0,t] 

and we take Yt — x[ n ^ ■ ■ ■ x\ X \ Assume the formula is true form n. The first term 
will be the sum over all partitions of {n + 1, n, . . . , 1} in which {n + 1} appears as 
a singleton, the second term will be the sum over all partitions of {n + 1, n, . . . , 1} 
in which n + 1 appears as an extra in some part of a partition of {n, . . . , 1}. In this 
way we arrive at the appropriate sum over ty n +i- 



Corollary 3. The inversion formula for off-diagonal integrals is 
I dX™ • • • dX™ = $>(.A) / dX% n) ■ ■ ■ dX% } 

J[o,tr J[o,t] N{v) 



where /j, (A) is the Mobius function introduced earlier. 

1.1. Multiple Martingale Integrals. We suppose that X t is a martingale 



so that, in the particular, E 



S[o,t) ^ ( s ) dX » 



where 4> is any adapted inte- 



grate function. In general, multiple integrals with respect to X will not have zero 
expectation, however, this will be the case for the off-diagonal integrals: 



E 



J[o,tr 



dX { t n) ---dX. 



(1) 



0. 



This property is in fact the main reason for introducing off-diagonal integrals in 
the first place. 

In the special case of classical processes we can employ the commutativity to 
write 



L 



dXi n) •••dAf ) = 



1 n\ 7[o, 



dX, 



(n) 



■dX. 



(1) 



'A"(t) n - J[0,t] r 

We define a family of random variables E Xt by 



n>0 



dX^ ] ■■■dX t 



[O.t] 



(1) 



for t > 0. We refer to process 1 1— > E * as an exponentiated martingale. 



1.2. Wiener Integrals. Let us take X t to be the Wiener process W s . We 
will have 



dW, 



(«) 



dW (n ] .■■■dW (1) (U 



Here we have dW s dW s = ds: this means that if we sum over partitions A in the 
inversion formula then we need only consider those consisting of singletons and 
pairs only. We will then have E (A) = n\ + 2n,2 = n and N (A) = n\ + n,2 = n — n2 ■ 
The combinatorial factor is F (A) = 1 . 
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Now there are ( ™ ) = ( 2 ™ ) ways to choose the singletons and ^"^Ji ways to 
choose the pairs. This yields 

- [«/2] 

- £(-■)— 2 .. n2 , ( :_ 2 „ 2) , w)- 2 "'" 

where _ff„ (x) = X)L"fo' ( — -^jA-—^ x n ^ k are the well-known Hermite poly- 
nomials. The implication that J A „ (t) dVK t ( ™ ) • • • dW^ = jjt n/2 H n is a result 
due originally to Ito . 

Using the relation Y^=o ^i,H n (x) — exp (xt — t 2 /2) for generating the Hermite 
polynomials, we see that the exponentiated random variable corresponding to the 
choice of the Wiener process is E zWt = e zWt ~ z */ 2 . 

More generally, if we take W t (/) = Jj t j / (s) dW s for / square-integrable then 
its exponentiated random variable is 



E Wt{f) = exp J W t (/) ~\f f M 2 ds 



1.3. Compensated Poisson Process Integrals. Let N t be the Poisson pro- 
cess and Y t = N t — t be the compensated process. We have the differential rule 
(dN t ) p — dN t for all positive integer powers p. The process Y t is a martingale 
and we have (dY t ) p = dN t — dY t + dt for all positive integer p with the obvious 
exception of p = 1 . 

It is convenient to replace sums over partitions with sums over occupation 
numbers. This time we find 

r- ^ E(n)=n n2+ 2„ 3+ 3„ 4+ ... ^ 

T [0 , tr dY ^'" dY ^^ ^ W (2)" 2 (3)" 3 • • • ni\n 2 \n 3 \ ■ ■ ■ ^ ^ + ^ ^ 
Again we use a Yl < — ¥ II 2 trick! 



/_^\n2+2n 3 +3ri4+--- 2 , r j 1 +2n 2 +37i 3 + 



k=l Kn k =Q K v 

= e -*J]exp (-lf-^^+t) 
fe=i ^ 

and using J2kLi (~^-) k 1 IT = m 0- + z ) wc cnc ^ U P w ith 

= e~ zt (1 + z) n+t . 

This means that, in terms of the Poisson process N, the exponentiated random vari- 
able associated with the compensated Poisson process Y t is E Yt (z) = e~ zt (1 + z) Nt 
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and we have 

/ dY^.-.dY^ =C n (N u t) 

where C n (x, i) are the Charlier polynomials determined by the generating relation 
J2 n >o frfC™ ( x > *) = e ~ Zt (! + Z ) X - Explicitly we have 

2. Ito-Fock Isomorphism 

Here we examine stochastic processes and exploit the fact that they can be 
viewed as commuting one-dimensional quantum fields (with base space typically 
taken as time). It should be remarked that the Wiener and Poisson process of clas- 
sical probability can be naturally viewed as combinations of creation, annihilation 
and conservation processes introduced by Hudson and Parthasarathy |10| as the 
basis of a quantum stochastic calculus which extends the Ito calculus to operator- 
valued processes. In fact the Ito correction can be alternatively understood as the 
additional term that arises from Wick ordering of stochastic integrals with respect 
to white noise creation/annihilation operators, see [5], for instance, and also [7|- 

2.1. Ito-Fock Isomorphism. Let X t be a classical martingale with canoni- 
cal probability space (flx,J~x,P) and let \)x — L 2 (fix, Tx, P)- We consider the 
function F (t) := E [Xf~\ and this defines a monotone increasing function. We shall 
understand dF (t) to be the Stieltjes integrator in the following. 

It turns out that our considerations so far allow us to construct a natural 
isomorphism between l)x and the Fock space r + (l? (R + , dF)), see e.g. |14j . 

For fe L 2 (K + , dF), we define the random variable X (f) := Jj Q ^^(s) dX s and 
X t (f):=X(l m f). 



Lemma 16. Let z\ n) (/) = / f . dX tn (/)■•■ X tl (/) th 



A»(t)' 



E 



Z[ n) {f) (g) 



1 


[L 







f(u)g(u) dF(u) 



Proof. For simplicity, we ignore the intensities. Let z[ n ^ = J 4 „, t , dXt n ■ ■ ■ dXt 1 

whenever n > 0. Next suppose that n 



then we have E [z t (n) zi 0) j = E \z[ n) 
and m are positive integers, then 

y(m) 



E 



zrz y 



= E 



dx u zi n ~v / dx v zi m ~^ 



dF[u) E 



7 (n-\) y(m-l) 



and we may re-iterate until we reduce at least one of the orders to zero. We then 
have 



E 



z (n) z (m) 



tAs 



dF{u n ) / dF (u„_i) • • • / dF( Ul ) 



= -F(tAs) n S ri 
n! 
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The proof with the intensities from L 2 (1R + , dF) included is then a straightforward 
generalization. | 



Theorem 11. The ffilbert spaces h x = L 2 (tlx, Fx, V) and T + (L 2 (R+,dF)) 
are naturally isomorphic. 

Proof. Consider the map into the exponential vectors (|1 .4|> given by 

E*W - e (/) 

for each fg L 2 (R + , dF). We know that the exponential vectors are dense in Fock 
space and in a similar way the exponential martingales E x ^ are dense in t)x- The 
map may then be extended to one between the two Hilbert spaces. 
Unitarity follows from the observation that 



E 



E x U) E x(g) 



fin \ fa <if 

3 J [0, oo) J ^ 



which is an immediate consequence of the previous lemma. 



The choice of the Wiener process is especially widely used. Here we have the 
identification 

L 2 (il w ,Fw,^) = T+ (L 2 (R+,dt)) 
which goes under the name of the Wiener-Ito-Segal isomorphism. This result is one 
of the corner stones of Hida's theory |9 of white noise analysis. The same Fock 
space occurs when we consider the compensated Poisson process also |14j . 
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